Wealth Distributions in Asset Exchange Models 
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FIG. 1: Illustration of the asset exchange model 



P. L. Krapivsky and S. Redner 1 
'Department of Physics, Boston University, Boston, MA 02215 

How do individuals accumulate wealth as they interact economically? We outline the consequences 
of a simple microscopic model in which repeated pairwise exchanges of assets between individuals 
build the wealth distribution of a population. This distribution is determined for generic exchange 
rules — transactions that involve a fixed amount or a fixed fraction of individual wealth, as well as 
random or greedy exchanges. In greedy multiplicative exchange, a continuously evolving power law 
wealth distribution arises, a feature that qualitatively mimics empirical observations. 
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PERSPECTIVE 

The economy is a complex interacting system that re- 
sponds to a multitude of influences and extends over a 
wide range of monetary scales. As experience with fi- 
nancial crises continues to demonstrate, understanding 
how an economy develops and how it is influenced by ex- 
ternalities remains poorly understood. Basic questions 
about what causes financial crises and how to deal with 
them continue to be hotly debated, with little sign that 
a fundamental understanding is emerging [H-Q. 

What can statistical physics contribute to this discus- 
sion? Not much, if the goal is to predict the economy 
next year. However, statistical physics possesses power- 
ful theoretical tools that have proven useful in describing 
specific financial phenomena, such as the Black-Scholcs 
options pricing formula Q ■ There are many parallels be- 
tween statistical physics and economic phenomena, and 
physics-based modeling has helped facilitate conceptual 
developments in finance and economics (see e.g., |5|47|). 

In classic economic theories, humans, or companies, 
are considered as rational actors that respond determin- 
istically to external conditions. More recently, stochas- 
tic tools have been applied to the economy, particularly 
to financial modeling. The stochastic approaches that 
are conventionally employed arc Brownian motion and its 
generalizations. In physics, a similar approach was fol- 
lowed to describe non-dctcrministic systems, where the 
interaction between a particle and its environment was 
mimicked by noise, while interactions between micro- 
scopic entities (such as Brownian particles) were ignored. 
This development (associated with physicists like Ein- 
stein, Langevin, and Stratonovich, and mathematicians 
like Kolmogorov, Feller, and Ito) led to increasingly so- 
phisticated stochastic processes 
that is still active. 



ll| , a research thread 



Over the last 40 years a new approach that combines 
the stochastic behavior of elemental entities as well as 
their mutual interactions has emerg ed (see e.g., 13 
While a few-particle interacting system is hopelessly com- 
plicated and beyond the reach of analytical techniques, a 
dramatic simplification arises for a many-particle system 



because we can often make statistical predictions about 
its fate. That a macroscopic interacting system is simpler 
than its few-element counterpart has several appellations 
- the law of large numbers, ergodicity, etc. — and it 
justifies the utility of statistical physics and probability 
theory in attempts to understand economic processes. 

In this short review, we present an interacting many- 
agent asset exchange model that can be quantitatively 
analyzed using statistical physics tools. An agent could 
be a single person or a self-contained economic entity, 
such as a company. In this model, the interaction be- 
tween two agents results in a redistribution of their as- 
sets. We regard an asset as any economic attribute - 
cash, goods, or other materials — that contributes to 
overall individual wealth. 

The macroeconomy is viewed as the result of a large 
number of asset exchanges between randomly-selected 
pairs of agents. Through these exchanges a global wealth 
distribution develops, and we want to understand how 
generic features of this distribution depend on the nature 
of the exchanges. The notion that the wealth distribu- 
tion is driven by two-person exchanges appears to have 
been first considered in the economics literature by An- 
gle 14 , lj| . In the physics community, this approach was 
introduced by Ispolatov et al. Hill, a nd related perspec- 



tives on this subject include Refs. [17H2 1| . A comprehen- 
sive review of this research topic is given in 1221 and an 
engaging non-technical exposition appears in [23l |. 



2 



ADDITIVE EXCHANGE 

As a preliminary, we first study additive exchange pro- 
cesses, in which a fixed amount of asset is exchanged be- 
tween two agents, independent of their wealth before a 
trade occurs. At the outset, we have to determine how to 
treat agents whose wealth reaches zero as a result of many 
unfavorable trades. We treat such penurious agents as 
economically "dead" , so that they no longer participate 
in the evolution of the wealth. Mathematically, this rule 
corresponds to imposing an absorbing boundary condi- 
tion on the density of agents of zero wealth. An alterna- 
tive is to impose a reflecting boundary condition at zero 
wealth, so that all agents continue to economically inter- 
act, even if they have no wealth [13]. In this latter case, 
the wealth follows the Boltzmann distribution of equilib- 
rium statistical mechanics, with an effective temperature 
equal to the average amount of wealth per agent. 

Under the condition that bankrupt agents are elimi- 
nated from further economic activity, we determine the 
consequences of: (i) fair transactions, where either agent 
is equally likely to profit in an interaction, and (ii) greedy 
transactions, in which the richer agent profits in an inter- 
action. For simplicity, each agent is assumed to possess 
an integer- valued amount of assets and that one unit of 
asset is transferred between traders in each interaction. 



Fair Transactions 

In a fair exchange, the wealth of two agents evolves 
as (j, k) — > (j ± 1, k =F 1); the direction of the exchange 
is independent of their starting wealth. The wealth dis- 
tribution evolves by selecting two agents at random who 
exchange one unit of wealth and repeating this elemental 
step ad infinitum. We assume that all agents are equally 
likely to interact with any other agent (corresponding to 
the mean-field limit in statistical physics). 

In this limit, the evolution of the wealth distribution 
is described by a master equation that accounts for the 
changes in wealth in each microscopic interaction be- 
tween agents. Let Cfc(i) be the density of agents with 
wealth k. In random additive exchange, the master equa- 
tion is 



dck at r 
~dT = N [Ck+1 



Cfc-i - 2c k ] 



(1) 



where N(t) = J2 k>1 Cfc(t) is the density of economically 
viable agents. The first two terms on the right-hand 
side account for the gain in c k due to the transactions 
(j, k + 1) — > (j + 1, k) and (j, k — 1) — > (j — 1, fc), respec- 
tively, while the last term accounts for the loss in c k due 
to the transactions (j, k) — > (j ± 1, k =p 1). Since these 
transactions require the presence of an agent of wealth k 
or k ± 1 and an agent of arbitrary wealth, all terms on 
the right-hand side involve N times a concentration. The 



density of agents with a single unit of wealth evolves by 
dci/dt = N(c2 — 2ci); this equation may also be written 
in the same form as Eq. (JXJ) by imposing the absorbing 
boundary condition Cq(T) = 0. 

Introducing the time-like variable, T = j Q dt'N(t'), we 
reduce Eq. J]]) to the discrete diffusion equation 



dck 
dT 



Ck+i + Cfc-i - 2c fc , 



(2) 



which may be solved for any initial condition [HI [l3| . 
When all agents start with unit wealth, c k (0) = 5k,i, we 
may account for the absorbing boundary condition by 
augmenting the initial condition with an "image" con- 
tribution due to agents with initial wealth —1; that is, 
Cfc(0) = 5 k ,i — <5fc,-i- The solution to Eq. © subject to 
these initial conditions is 11311 



c k (T) 



-2T 



[h-i (2T) 



(3) 



4+1 (2T)] . 

where /„ is the modified Bcsscl function of order n. Cor- 
respondingly, the total density of active agents N(T) is 

N(T) = e- 2T [J (2T) + h{2T)\ . (4) 

the asymptotic behaviors of 



In the limit T - 
Eqs. © and (gj) arc 

k 



k 2 /4T 



N ~ (ttT) 



-1/2 



(5) 



VIttT 3 

These asymptotics apply, up to an overall factor, to all 
initial conditions that decay sufficiently rapidly with k. 
An important feature of §5§ is the emergence of scal- 
ing: the distribution c/ £ (T) depends on the scaled wealth, 
k/VT, rather than separately on the variables k and 
T. Similar scaling behavior arises in numerous inter- 
acting particle systems fl3| . Normally, scaling is pos- 
tulated and then verified analytically or numerically. For 
asset exchange, we deduce the validity of scaling from 
the exact solution. We now express the asymptotic 
solution ([5]) in terms of the physical time t by using 



t(T) = £ dT'/N(T') 



Ck 



— exp 
3t 



2 V 7rT 3 to eliminate T and give 

1/3 



3 

/ 7T x 1 / 3 k" 



V144V 



t 2/3 



TV 



2 



(6) 

The number of viable agents decreases as t and their 
typical wealth grows as i 1 / 3 . While this model is not 
realistic, it illustrates the efficacy of a statistical physics 
perspective in solving an interacting many-body system. 

Instead of removing bankrupt agents, let us provide 
each of them with 'welfare' of a single unit of asset. In 
this case, the economically viable population density is 
always N = 1 and the master equation for the wealth 
distribution simplifies to 

dck 

— = c fe+ i + c fc _i - 2c k fc>2, 

dt (?) 

dci 

— = c 2 - ci , fc = 1 . 
dt 
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We can extend the first of these equations to all k and 
also subsume the equation for c\ by choosing the initial 
condition Ci_fe(0) = c k (0), with ci(0) = c (0) = 1, and 
c fe (0) = for fc ^ 0, 1. The solution to © subject to this 
initial condition is 



Ck 



[Jfc_i(2t) + I fc (2t)] 



(8) 



Because of this injection of assets to destitute agents, the 
total wealth density of the population, M = J2k>i ^ c fc> 
grows with time as 



M = e 



-21 



^fc[4_x(2t)+4(2t)] 



k>l 



as t — > oo. In this toy model, the rate of welfare expen- 
diture to keep everyone solvent decreases with time! 

Greedy Transactions 



where C = dC/dx. The scaling function must satisfy 



dx C{x) = 1, and 



dx x C(x) = 1, (12) 



that follow from N — c k(t) — N J dxC(x) and setting 
the (conserved) wealth density to one, J^k ^ c k(t) = 1- 

Equation (fTTj) is soluble by elementary techniques [l6| , 
and the asymptotic wealth is simply the step function 



c fc (t) 



1/(2*), k<2y/t, 
0, k > 2\/t, 



(13) 



while the density of active agents decays as N(t) = t~ x / 2 . 
In greedy exchange, the number of viable agents decays 
faster than in random exchange and the population is 
slightly wealthier, with the average wealth growing as 
t 1 / 2 rather than as t 1/3 . 



In greedy exchange, the richer agent is exploitative and 
always takes one unit of wealth from the poorer agent in 
each interaction, as represented by (j, k) — > (J + 1, k — 1) 
for j > k. The densities Cfc(i) now evolve according to 



dek 
dt 



fe-i 



Cfc-i > Cj + c k+ i 



2^ c j 



c fc (iV + c fc ). (9) 



The first term on the right accounts for the gain in Ck 
due to an agent with wealth k — 1 taking one wealth 
unit from a poorer trading partner. Similarly, the second 
term accounts for an agent with wealth k + 1 losing one 
unit of wealth to a richer trading partner. The last term 
accounts for the loss of Ck when an agent of wealth k 
trades with anyone; the extra factor of c& accounts for 
the loss of both agents of wealth k when two such agents 
interact. 

While this set of non-linear equations appears in- 
tractable by exact methods, they are readily amenable 
to a scaling analysis 13J . We first re- write Eq. ([9]) as 

dek 

— = -c k {c k + Cfc+i) + N(c k -i - c k ) 

« (10) 

+ (Cfe+l -Cfc_i)2^Cj, 

and make the scaling ansatz Ck(t) ~ k~ 2 C(k/k t ), where 
fc*(t) is the typical wealth of each agent. That is, the 
wealth distribution at different times is invariant when 
wealth is measured in units of the time-dependent typical 
wealth. The prefactor k~ 2 ensures that the total wealth 
of the population, kck(t), is conserved, while the con- 
dition J2k c k(t) = N{t) gives K(t) - l/N(t). Substitut- 
ing now the scaling form Cfe(i) ~ N 2 C(x), with x = kN, 
in Eq. (|10p and taking the continuum limit gives 



C(0)[2C + xC'} =2C 2 +C 



1 - 2 



dyC(y) , (11) 



MULTIPLICATIVE EXCHANGE 

While additive exchange provides instructive warm- 
up examples, multiplicative exchanges, where a fixed 
fraction of the current wealth of one of the agents is 
traded, are economically more realistic. For example, 
investment returns are generally quoted as percentages 
rather than absolute amounts. A trade now has the form 
(x, y) — > (x — ax,y + ax), with < a < 1 the fraction 
of the loser's assets that are gained by the winner. By 
multiplicative exchanges agents can never go bankrupt, 
but they can become arbitrarily poor. 



Fair Transactions 

If an agent gains or loses with equal probabilities in a 
transaction, the wealth distribution evolves as 

(Jcix) 

—jjf - = \H d v dz c (y) c ( z ) x [ — S(x — z) — S(x — y) 

+5(y(l - a) - x) + 5(z + ay - x)] . (14) 

The delta functions cleanly indicate the origin of the 
various terms in this equation. For example, the first 
two terms on the right account for the loss of agents 
of wealth x due to trades with any other agents. The 
next two terms account, respectively, for the gain in 
c(x) due to the exchanges (jz—,y) — > (x, y + jz^) and 
(y, x — ay) — > (y(l — a),x). Integrating over the delta 
functions, the master equation becomes 



dc(x) 
dt 



< x ) + 3(i=sy c (r^) 



+ ^ [ X dyc(y)c(^), 
Jo 

where we set the (conserved) total density to one. 



(15) 
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Equation (jT5|) is daunting, and it is simpler to study 
the evolution of the moments, M n {t) = L dxx n c(x,t), 
that quantify the wealth of a typical agent. It is straight- 
forward to verify that the first two moments, the popu- 
lation Mo and the wealth density Mi, are conserved; we 
choose AIq = 1 and Mi = M without loss of generality. 
More interesting behavior arises for the second moment 
equation 



dM 2 (t) 
dt 

whose solution is 
M 2 



-a(l - a)M 2 (t) + aM 2 



M 2 (t) = 



1 - a 



M 2 (0) - 



M 2 



1 — a 



,-a{l-a)t 



(16) 



All moments beyond the second also converge to non-zero 
steady-state values. This steady state arises because the 
wealth of a rich agent substantially diminishes in a los- 
ing multiplicative exchange, but its wealth increases only 
slightly in a winning exchange. Conversely, a poor agent 
suffers a slight loss in a losing exchange but can gain sub- 
stantially in a winning exchange. These two countering 
outcomes tends to move all agents toward a middle class. 



Greedy Transactions 

When only the richer agent gains in an exchange, the 
master equation is now 



<9c(a 



dt 



POO 

-<*=)+ / 

Jx/(l-a 



dy c{y) 



/(!-«) 

x/(l+a) 



(17) 



Numerically, we find that the resulting wealth distribu- 
tion is a power law (Fig. [2]), with most of the population 
impoverished. Pervasive impoverishment arises because 
greedy exchange causes the poor to become poorer and 
the rich to become richer, but wealth conservation forces 
there to be many more poor than rich agents. In the long- 
time limit, a small fraction of the population possesses 
most of the wealth. 

An exact formal solution to Eq. (|17j) is [l6| 



c(x, t) 



A 

xt 



with A 



1 



ln(l - a) 



(18) 



This distribution is pathological, however, because posi- 
tive the moments M n (t) of this distribution are divergent. 
Thus Eq. (|18j) can only apply within an intermediate scal- 
ing regime xi(t) < x < X2(t), a restriction that leads to 
finitcness of all the moments. To determine this scaling 
region, we use Eq. (|18p to compute the moments and 




FIG. 2: The power-law wealth distribution c(x) of greedy 
multiplicative exchange with a — 0.5 on a double logarithmic 
scale for times t — 1.5™, with n = 7, 10, 13, and 16. 



obtain: 



M (t) 
Mi(t) 



X2 dxc(x,t)~^M 
Ax 2 



dx x c(x, t) 



Xl 



t 



(19) 



Since Mq = 1 and Mi arc constants, wc infer that 
xi(t) — e~ l / A = (1 - a) 1 and x 2 {t) oc t. These cut- 
offs correspond to the wealth of the poorest and richest 
agent, respectively, in the population. It is only within 
these ranges that the wealth distribution is a power law, 
as shown in Fig. [2] 



DISCUSSION 

Asset exchange represents a parsimonious mechanism 
for the gain and loss of individual wealth in an econom- 
ically active population. In spite of the obvious short- 
comings of considering only this single factor among the 
myriad of influences on individual wealth, asset exchange 
models lead to a rich array of wealth distributions. For 
additive asset exchange, the wealth distribution can be 
explicitly derived for a variety of microscopic exchange 
rules. For greedy multiplicative exchange, where the 
richer agent always gains in an interaction, a scaling- 
based approach indicates that the wealth distribution has 
an evolving power law form, c(x, t) cx l/{xt). 

Power-law distributions occur in the high-wealth tail 
of the wealth distribution in various economies, with 
the associated exponent in the range of 1.6-2.2 (see 



Refs. [21j, |24j). As alluded to in the introduction, 
a variety of stochastic models, where agent undergoes 
an independent stochastic process, have also been in- 
voked to argue for this power law (25rf29j . In contrast, 
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greedy multiplicative exchange is based on a combina- 
tion of stochasticity and microscopic interactions be- 
tween agents. There are many directions in which as- 
set exchange models have been extended to make them 
more realistic; recent work along these directions can be 
found in Ref. (30l436j . Specific examples of such addi- 
tional elements include the incorporation of the saving 
of assets 1^, 21, 33], speculative trading [38| , and other 
forms of wealth redistribution. The notion of exchange 
of assets has also been applied to construct a migration 
model for the distribution of city sizes |3{|. It should 
prove interesting to examine the role of such redistribu- 
tion mechanisms in the ideologically-free setting of sta- 
tistical physics modeling. The underlying assumption of 
conserved assets in an exchange neglects the possibility 
of wealth growth because of the exploitation of a natural 
resource, technological developments, or by both agents 
benefiting in exchanges. These are issues that appear 
ripe for further development. 

We thank Slava Ispolatov for his initial collaboration 
on this project. 
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